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Abstract 

A variational method is used to analyse compact (7(1) gauge theory in 2 -f- 1-dimensions 
at finite temperature, T, weak coupling, g and where the fundamental magnetic monopoles 
have magnetic charge 2nn/g. The theory undergoes a critical transition from a confining 
phase at temperatures below T c = 2g 2 /n 2 7r to a deconfined phase at temperatures above T c . 
The free energy and all its derivatives are continuous at T c , indicative of the BKT phase 
transition. The relevant gauge-invariant correlation functions decay exponentially at large 
distances. The spatial Wilson loop obeys the area law at all finite temperatures, even for 
the non-compact theory. The case n = 2 corresponds to the compact (7(1) theory considered 
as a low energy effective theory for the spontaneously broken Georgi-Glashow model. The 
results in this case agree with those derived previously for compact (7(1) in this model using 
dimensional reduction of the Lagrangian. 

PACS Numbers: 12.20.Ds, ll.10.Kk, ll.10.Wx, 11.15.Tk 
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1 Introduction 

A recent paper jj]] uses a variational method to study the deconfinement transition in QCD at 
finite temperature. The method mimics the Rayleigh-Ritz variational method in the Schrodinger 
formulation of quantum mechanics. There the standard procedure is to form an ansatz for the 
ground state wavefunction, parameterized by some free parameters, and to minimise the expec¬ 
tation value of the Hamiltonian with respect to the parameters. This provides an upper bound 
for the ground state (vacuum) energy. The method at finite temperature is analogous: in the 
canonical ensemble formulation of quantum statistical mechanics, one forms an ansatz for the 
density matrix, with free parameters, and minimises the expectation value of the Helmholtz free 
energy. This provides an upper bound for the free energy at a given temperature. 

This paper carries on the theme of Q, extending the variational analysis to the compact (7(1) 
gauge theory in 2 + 1-dimensions at finite temperature. This theory was shown to exhibit confine¬ 
ment at zero temperature, even at weak coupling, by Polyakov many years ago Q . Essentially, the 
confinement is due to gauge field vortex defects which act in the partition function as a gas of inter¬ 
acting magnetic monopoles, causing Debye screening of the gauge field, a non-perturbatively small 
mass gap for the photon, an area law for the Wilson loop and confinement of electric charges. This 
was confirmed by a variational analysis More recently, the theory has been studied at finite 

temperature as the low energy limit of the Georgi-Glashow model (5(7(2) gauge field cou¬ 
pled to adjoint Higgs). At weak coupling g and low temperature T, the theory can be understood 
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as a two dimensional Coulomb gas, undergoing the BKT phase transition ]7j, at temperature 
= i? •0 However, because the interaction strength in the Coulomb gas picture is temperature 
dependent, the phase structure is reversed , with a dilute plasma of unbound monopoles below T c 
and a gas of bound magnetic dipoles above. This further implies the deconfinement of electric 
charges above T c . As we shall see below, this phase transition is indeed apparent in the variational 
analysis at finite temperature: for the compact 17(1) theory considered on its own, with fundamen¬ 
tal magnetic monopoles j of magnetic charge 27 Tn/g, there is a phase transition at temperature 
T c = 2g 2 /n 2 n (eq’n |67| ) from a confined phase (with finite photon mass gap) to a deconfined 
phase (zero mass gap). The mass gap vanishes continuously, indicative of a critical transition. 
Furthermore, one can show that the free energy and all its derivatives are continuous at the phase 
transition (^9|). This is strongly indicative of the BKT phase transition. For the compact 17(1) 
theory considered as the low energy limit of the spontaneously broken Georgi-Glashow model, in 
which the ’t Hooft-Polyakov monopoles have fundamental charge 47 r/g 0. the phase transition 


is at T c = as expected. 

The analysis proceeds as follows. Since a variational method is only as good as its ansatz , it 
is prudent to begin by studying the form of the density matrix in simpler theories, starting with 
the simplest of all, the harmonic oscillator. It is trivial to construct the density matrix in the 
canonical ensemble in the energy basis, less so in the co-ordinate basis. But this is precisely the 
basis one needs for applying the variational method to theories where the spectrum is unknown. 
In fjl], it is shown that the density matrix for the harmonic oscillator takes the form of a Gaussian. 
This generalises to free scalar field theories upon replacing co-ordinates by fields. 

Section || begins by outlining the details of the variational method. The variational ansatz 
is taken to be a Gaussian, with arbitrary kernels. The minimisation of the free energy is then 
discussed; it is shown that the minimisation for the harmonic oscillator picks out the density 
matrix constructed previously in §|j. 

Section discusses compact C/(l), dealing first with the requirement of gauge invariance. A 
straightforward modification of the ansatz for the density matrix (following Q]) projects onto 
the gauge-invariant sector of states in the theory. Thus armed, we proceed to the variational 
analysis of compact U(l), where things are complicated by the singular gauge defects acting as 
magnetic monopoles. However, at weak coupling, the monopole fugacity (and photon mass) are 
non-perturbatively small and one expects that the density matrix will not be significantly changed 
from the non-compact case. This is indeed what happens: upon minimisation of the free energy 
one finds a solution for the kernels (|5^) with a photon mass-squared proportional to the (non- 
perturbatively small) monopole fugacity (58). It is further shown that (as the coupling tends to 
zero) the variational method reproduces the expected behaviour of non-compact 17(1), which is 
after all just a gas of free photons with a single massless transverse degree of freedom (in 2 + 1 
dimensions). 

Once the variational solution for the kernels has been obtained, it is straightforward to analyse 
the phase structure of the theory (§^). Substitution of the kernels into the definition (33) of 
the monopole fugacity z yields an implicit equation for the fugacity (|63|). This equation has 
two solutions: the first is the trivial one z = 0; the second has z non-perturbatively small for 
temperatures below T c , zero at T c and non-perturbatively large above T c . The free energy is 
minimised by the latter solution below T c and by the former solution above T c . Finally, it is 
shown that all derivatives of the fugacity (and hence the free energy) are continuous at T c , a 
signature of the BKT phase transition. 

In section |^, the compact gauge invariant correlation functions and the spatial Wilson loop 
are calculated. The Wilson loop has the area law at all non-zero temperatures, even for the 
non-compact theory. This is confirmed in the Lagrangian approach to finite-temperature field 
theory. 


1 A more recent work 'j including the dynamics of W 1 Bosons in the G—G model has shown that there is a 
phase transition at an even lower temperature corresponding to the unbinding of IV’ 1 . This does not concern us 
here: we are interested only in the compact 17(1) aspect. 

2 By fundamental we mean that the smallest magnetic charge is 2nn/g-, integer multiples of this are of course 
allowed. 
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2 The harmonic oscillator revisited 

Construction of the density matrix for the harmonic oscillator in the energy basis is trivial. Let the 
normalised energy eigenstates of the Hamiltonian Tt be | n) with energy E n = (n+^)u> (henceforth, 
h = m = ks = 1)- In the canonical ensemble, the density matrix (normalised to the partition 
function) is then 

p = e~ H,T = Y, n e~ En/T \n)(n\. (1) 

In the co-ordinate basis, this becomes 

p(x, x') = Yi n e~ En t T {x\n){n\x'). (2) 


The Schrodinger equation is 


with solution 


1 d 2 l 2 2 f 1 , 
W + 2 W:C -{ n+ 2 |W 


{x\n) = 0, 


(3) 


H n (u > 2 x)e , 


^7r(2 n n!) 2 / 

where H n is a Hermite polynomial. Substituting in and using the formula [ ftl| 

1 / Z\ U „ , , „ , , 1 Z r. 9 9 . -1 

[(x + V )z - 2 xy \, 


/ z\ n 

H n (x)H n (y) 


1 ^ 
- ex P~I- 2 

1 — Z A 


VI - z 1 


(4) 


(5) 


one obtains 


/ (u\h e “/ 2T u 

p( x,x)=[-) ^====exp 


2 e 

7 TJ V! - e- 2 “/ T ^ 2(1 - e ~ 2u / T ) L 


(1 + e~ 2uj/T )(x 2 + x 2 ) - 4e““ /T W]. (6) 


As a check, taking the trace one obtains 


trp = dx p(x , x) = 


o/2 T 


1 — e - “/ T ’ 


(7) 


the standard partition function as required. It is clear that this result generalises to free scalar 
fields, which consist of many such modes labelled by momentum p with dispersion relation w(p). 
The density matrix can be written as a product over all the modes, or equivalently as a sum 
over modes in the exponent. Swapping co-ordinates x , x for scalar fields A , A and dropping the 
normalisation factor, we can write the density matrix for scalar fields as 


p[A,A] 



dp u> 

2n (1 — e~ 2u, /T) 


(1 +e- 2 “ /T )(A 2 +A' 2 ) 


4 e-“ /T AA' 


( 8 ) 


We note two points at this stage. Firstly, at zero temperature, the density matrix reduces to 




(9) 


in which terms coupling A to A! vanish. But this must correspond to a pure state: the ground 
state. So, as a general rule, terms coupling A to A! correspond to finite temperature thermal 
disordering effects. 

Secondly, we see that the density matrix for free theories (Le. the harmonic oscillator) takes a 
Gaussian form.^ Moreover, it is possible to establish the converse: all Gaussian density matrices 

3 This is fortuitous, since then the ensemble average of observable O, which is trpO/trp, is a tractable Gaussian 
functional integral. 
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correspond to free theories. To do this write the most general Gaussian density matrix (which 
must be Hermitian and therefore symmetric in the real case) as 


p(x,x ) = exp— — [xG x x + xG x x — 2xHx]. 


Then reverse the arguments of equations (0]3), substituting ( G,H ) for (w,T) via 


( 10 ) 


,_i 1 + e" 2 “/ T 


G" 1 =! 


H = 2u 


\ _ g-2oj/T 

e“/ T 


| _ e —2 {jj/T ' 


( 11 ) 


It is thus evident that the most general Gaussian density matrix © can be formally written as 
the operator 


P = 


o-'H/T 


( 12 ) 


where Ti is the Hamiltonian for an harmonic oscillator, with (u>,T) some functions of ( G,H ) 
obtained by inverting (©• Note that this argument holds independent of thermal equilibrium 
considerations: (w, T) are to be regarded merely as a re-parameterisation (actually a diagonalisa- 
tion) of ( G,H ). This one-to-one identification between Gaussian density matrices and harmonic 
oscillators will turn out to be crucial in determining the entropy corresponding to (|lo|). 


3 The variational approach 

In the preceeding section, the density matrix has been derived for the harmonic oscillator and for 
the generalisation to free scalar fields. Now, for the purposes of pedagogy, let us derive the same 
results using the variational method. A suitable ansatz (the most general tractable one) is the 
Gaussian form 

p[A,A] = exp -~[AG- 1 A + A'G- 1 A' - 2AHA'}, (13) 

where we employ a matrix notation for integrals 

AGHA= ( dxdydz A{x)G{x — y)H(y — z)A(z) (14) 


etc. and the functions G and H are arbitrary kernels with respect to which the free energy will be 
minimised. 

The Helmholtz free energy is F = U — TS where U is the energy and S is the entropy. 
Remembering that p is no longer normalised, one has 


U = 


tvpTi 
tr p 


and 


s = tr[plog(p/trp)] = _trplogp +logtrp 


trp 


trp 


Here, H is the Hamiltonian 


H[A,A] = -5(A-A) 


S 2 

6 A 2 


+ Auj 2 A 


(15) 


(16) 


(17) 
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The ensemble average of the energy ([Is]) with respect to (jTo|) is a simple Gaussian integral. One 
obtains 


U = 


2tr p 
1 

2tr p 
1 


DAD A S(A-A) 


6 2 

6 A 2 


+ Auj 2 A 


p[A , A'] 


J DA [trG -1 - AiG- 1 - H) 2 A + Auj 2 A\ p [A, A] 


= ^tr(G _1 + H + w 2 (G _1 - H)- 1 ). 


In momentum space 


u =\ / ^ [G- 1 + H + w 2 (G -1 - JT)" 1 ] , 


(18) 


(19) 


where V is the spatial volume. How to evaluate the entropy © in general is less clear. Indeed it 
is not apparent how to define the logarithm of p. However, for a Gaussian density matrix we can 
use the harmonic oscillator correspondence proven in the previous section. In terms of (ui, T), the 
entropy is given by 


d 


cue 


—uj/T 


S = af TloeZ= ^l ° g ( 1 - e ~“' T ) + T(1 _ e -w/ T y 

which can be re-written in terms of G and H (in fact in terms of just h = GH/2) via (© as 
S — log 


( 20 ) 


2 h 

1 

' 1 - ( 1 - 4 / 1 2 ) 1 / 2 ' 

( l-(l-4 h 2 ) 1 ' 2 \ 

(1 — 4/1 2 ) 1 / 2 — (1 — 2h) 

log 

2 h 

i ((1 - 4*, 2 ) 1 / 2 — (1 - 2*,)) J 


( 21 ) 


The free energy in terms of G and h is, from (|l^) and (pl|), 


F = - [G _1 (l + 2h) + w 2 G(1 - 2k)- 1 } 


f ^ 

2 h 


1 — (1 — 4ft, 2 ) 1 / 2 

/ 1 - (1 - 4h 2 ) 1/2 \\ 

1 

fil-Ah 2 ) 1 / 2 - (1 — 2/l) 

log 

2 h 

\ (1 — 4/j 2 ) 1 / 2 — (1 — 2h)) ) 


Minimising with respect to G, one obtains 

G -1 = u>{l -Ah 2 )- 1/2 . 
Minimising with respect to h , one obtains 

_, 2 T 


0 = G _i + w 2 G(l — 2h)~ z + 


■log 


1 - (1 -Ah 2 ) 1 ' 2 


2 h 


(1 + 2/i) 1 / 2 (l — 2/i) 3 / 2 
Substituting the former expression for G into the latter expression yields the solutions 

1 + e~ 2uj/T 


( 22 ) 


(23) 


(24) 


G" 1 =u 


H = 2uj 


\ _g-2 u/T 

o-u/T 


\ _ g-2 u/T ’ 


(25) 


which are the kernels derived in section g(0 
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4 Compact U(l) gauge theory 


Gauge theories present an added complication: it must be guaranteed that gauge invariance 
is respected throughout. So physical observables must be gauge invariant and furthermore we 
must only consider states that belong to the gauge-invariant sector of the Hilbert space. For the 
variational method at T = 0 this is achieved by integrating over the gauge group in the ansatz 
for the wavefunction to project out gauge-invariant states |l^]. Correspondingly, in (j|) we should 
restrict the sum to gauge-invariant states, so we project the ansatz onto the gauge-invariant sector. 
Because the density matrix is bi-linear in the states, this requires two gauge projections. 

Consider a theory of a gauge field A M with gauge symmetry under a gauge transformation p 
such that A M —> A Working in the Hamiltonian formalism of the field theory with Aq = 0, the 
required gauge invariance can thus be achieved by modifying the density matrix ansatz (13) to 


p [A, A] = 


DpDp' exp — — 


Af G~ x Af + A'/G-^A/ - 2AfHA 


(26) 


Here i labels the spatial components of the gauge field. But if operator 0[A , A] is gauge invariant, 
one of the gauge group integrals in tr pO is redundant. It suffices to take 


p[A A } 


Dtp exp — 


1 ' 
2 . 


AtG~ At + A?G A? - 2 A Z HA* 


(27) 


Consider now applying the variational method to non-compact U( 1) in 2 + 1 dimensions. The 
relevant Z7(l) gauge transformation is Ai —> Af = Ai — dip/g , where g is the coupling and p 
is a regular function. We could now follow the variational method recipe given above verbatim: 
determine the energy with respect to the density matrix ansatz (|27| ) and form the free energy F, 
including the entropy of a single massless degree of freedom^] Then minimise F with respect to 
the variational kernels G and h. We choose not to do so here, because the result can be recovered 
from the compact theory discussed below by letting the coupling go to zero (in which case the 
magnetic monopoles decouple from the photons). One finds that the kernels are the same as for 
the single free massless scalar field, as one expects. 

The foregoing discussion has given us a detailed understanding of the density matrix and finite 
temperature variational method for what are essentially trivial free theories. It is now time to 
consider a less trivial theory: the compact 17 ( 1 ) gauge theory. In the Hamiltonian formalism, the 
compact theory differs from the non-compact theory in that it admits singular gauge transforma¬ 
tions p{x) which create vortex-like gauge field configurations which act in the partition function as 
magnetic monopoles. At zero temperature these monopoles form a dilute plasma, causing Debye 
screening of the 1 /( 1 ) field. Thus the photon has a mass gap, the Wilson loop exhibits an area law 
and the theory is confining. 

An important point to note is that the compact U( 1) theory is in a sense very similar at 
weak coupling to the non-compact U( 1 ) theory: the monopole plasma is extremely dilute, and the 
photon mass is non-perturbatively small. Thus one expects that the Gaussian variational ansatz 
( f27| ) will be a suitable one. 

The variational method at finite temperature followed below very closely parallels the analysis 
at zero temperature given in to which the reader should refer for details. The principle is 
essentially the same, but the presence of the additional kernel H and entropic effects complicate 
the algebra. 

We start with a short discussion of thermal ensemble averages. Consider first the simplest of 

4 Entropically speaking, free photons in 2 + 1 dimensions with a single transverse degree of freedom are equivalent 
to a single massless scalar degree of freedom. 
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all viz. tr p. This is of course just the normalization. Now 


trp = J DipDAi exp ~^[A i G~ 1 A l 
= J D^pDAi exp — A^' 2 {G~ 1 
= J DipDAi exp—tT. 


AfG^Af - 2A t HAf] 

H) A*' 2 + ^dMG- 1 + H)d t p 


( 28 ) 


As we shall see below, ensemble averages translate into the language of field theory as correlation 
functions evaluated with respect to Euclidean action W. 

To continue the evaluation of (p8|), split the gauge transformation ip into parts ip and <p v , 
corresponding respectively to regular (non-compact) gauge transformations and to a distribution 
of n_|_ vortices at points x a and n_ anti-vortices at points xp. Explicitly, 


<Pv = Z2=i n 6( x ~ x a) - £p=i n Q(x - x p), (29) 

where 6 is the plane polar angle and n is the charge (in units of 2n/g) on a fundamental monopole. 
Then the normalisation integral ( p8| ) can be evaluated as follows. The integrals over At and ip are 
Gaussian. Further, substituting the general vortex distribution for ip v we can write the ip v integral 
in standard fashion as the partition function Z v for a Coulomb gas. In all, 


trp = det [7 t(G 1 


H)~ l ] det 


4t rg 2 

d 2 


n i/2 


(G-'+H)- 


Zvi 


(30) 


where 


y _ yiCxD _ 

17 “ n +> n - =0 n + ln-\ 


n ati n /3=l / d X ad Xf} Z 


exp - 


4 g 2 


n_|_+n_ 


1 + H] + Y^p^D^xp - Xp>,G +H]-'£ a ,pD(x a -xp,G + H] 


(31) 


In this last expression for the vortex partition function, the vortex-vortex interaction is 

d 2 p G-' + H_ inx 


D(x, G 1 + H] = 8n 2 n 2 f 

Jo 


(27T) 2 


j.r 


the vortex fugacity is 


z = A 2 exp 


1 

"8g 


2 D(0 1 G~ 1 + H} 


(32) 


(33) 


and A is the UV cutoff in the theory. 

Note that, since the minimal kernels G and H will in general have temperature dependence, 
the monopole interaction D and fugacity z will also acquire T dependence. This is in accord with 
the dimensional reduction arguments of ||. 

After this brief diversion, let us get back to the discussion of thermal ensemble averages. 
Ensemble averages translate into field-theoretic language as correlation functions with respect to 
the ‘action’ W in (p8|). This action has three fields: A,, (p and ip v . The former two are free fields. 


Using the matrix notation of (14) and ignoring normalization factors (these cancel in practice), 


5 The notation D(x,G 1 + / 1 ] indicates that D is a function of the separation x and a functional of the kernel 
G -1 + H. 
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one can immediately write down the following correlators: 


(Ai) a — 0 dup, 
2 g 


(AiA^A = (G- 1 - H)- 1 + -^diipdup, 

4 g z 


{dw)<p = 0 , 

(di<pdi<p)e =2g 2 (G~ 1 + H)- 1 . 


(34) 

(35) 

(36) 

(37) 


The vortex field correlators are evaluated as follows. We write them as correlation functions of 
the vortex density 

p(x) = T, a pn[6(x — x a ) — S(x — x@)}. (38) 

The two-point vortex density correlator is identical to that in ||, if one takes account of the 
modified form of z in (33) . The result is 

(pp) p = 2z + 0(z 2 ). (39) 


These results are all that is needed to determine the energy with respect to the ansatz 
The Hamiltonian for compact 17(1) is 

H = ±[E? + b 2 ]. 

As per the Hamiltonian for non-compact U{ 1), the electric contribution to the energy is 

However, the magnetic part is not the usual 

±B 2 = ±6(A-A')(e ij d i A j ) 2 . 

This is because B is not invariant under compact gauge transformations, which in general create 
sources of magnetic flux, i.e. the monopoles. We require a singlet part of B 2 which is invariant 
under compact gauge transformations and which we denote by b 2 . 

The normalised electric contribution to the energy in the thermal ensemble is thus 


(40) 


(41) 


(42) 


trp ^ 2 


1 


« 2 


^r = -2^J DADA '^ A - A '^ A ’ A ^ < 43 > 

Doing the functional differentiation and performing the Ai, (p and tp v integrations in turn, this is 
1 


2tr p , 

2trp , 

2trp , 
1 


DADp S(A - A') S^ti-G- 1 - (A.G" 1 - A^H^G^Ai - HA?) 


„-w 


DADip 
DADipDp 


tr(G _i + H) - —jdiip(G~ + HYdw 
4 5 


—W 


itr(G _1 + H) - + B) 2 dip v 


-w 


1 


tr(G _1 + H) - —(pDiiG- 1 + H) 2 ]p), 

where Z)[(G _1 + H) 2 ] is the vortex interaction for kernel (G _1 + H) 2 , cf. 
space, using (pp) p = 2z, one gets 


fa 'P\ E 1 _y 

tip 


d 2 p 

WT 2 


ha-' + H)-’- 

4 g 2 p 2 


(G-'+Hf 


(44) 

In momentum 

(45) 





















This has the same form as in the T = 0 calculation with the replacement G _1 —> G _1 + H\ the 
difference is that G, H and z will acquire T dependence upon minimisation at non-zero tempera¬ 
ture. 

Next consider the magnetic contribution to the energy. The compact 17(1) Hamiltonian con¬ 
tains 6 2 , a part of B 2 which is a singlet under compact gauge transformations. The existing 
literature on the Hamiltonian treatment of compact 17(1) atT = 0 suggests two ways in which to 
proceed here. The first || is to note that the ensemble averages of trpb 2 = tr pB 2 should coincide 
when taken with respect to the properly compact gauge-invariant density matrix, so it suffices to 
calculate the latter. There is a caveat however: because B 2 is not compact gauge invariant, we 
must use (^ 6 |) rather than (27) as the ansatz for the density matrix, since the latter form is only 
valid for operators which are fully gauge-invariant. 

The second method is to explicitly choose a singlet form for b 2 . In |Q|, the authors chose the 
‘mixed action’ 


b 2 = 



a cos gB + 


(!-«) 

4 


cos 2 gB , 


(46) 


which reduces at weak coupling to the standard non-compact magnetic term. At zero temperature, 
this theory contains a surprise: for a > 7 r 2 /4 there is no mass gap in the best variational state. So 
the different methods give different results at zero-temperature, presumably because the singlet 
parts are different. However, we note (but do not show here) that at finite temperature, the 
analysis of the theory with mixed action parallels the first method above: if there is a mass gap 
at T = 0 (a < 7 r 2 / 4 ), it will disappear at T c ; if there is no mass gap at T = 0, there is no mass 
gap anywhere. 

Following the dogma of the first method, one has that 


trp^- trp 
trp tr p 

1 


2 trp 


7 / DAD A Dip'Dip" B 2 exp - 


A? G _1 Af + A * G~ l A? - 2 A? HA l 


(47) 


// / 


Now make the change of variables A, A —> A v , A v and let ip = ip — tp , p = ip + ip such 
that 


trp^r- 1 


trp 2 trp‘ 


7 / DADipDr] 


di ( A o - ^ 


1 2 


- W 


(48) 


The term linear in 77 vanishes, and the quadratic term is independent of the variational kernels; 
we omit it. Then 


trp^ 1 

trp 2 trp' 

Making the change of variables A 


DADipDr] 

■ A~' p / 2 yields 


e?/ c)j 





-w 


(49) 


^ = J*r (50) 

Altogether, the ensemble average of the internal energy, the sum of electric ([fo]) and magnetic ( |5C| ) 
parts, is 


U = V 


cPp 
4(2t r ) 2 


G 


-1 


H - An2x l Z (G- 1 + H) 2 + p 2 {G~ 1 ^ H)- 1 
9 P 


(51) 
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Again it is unclear how to determine the entropy from first principles]^ We conjecture that the 
entropy written in terms of the kernels is the same as for a free theory, as in ( pT| ) . Why should 
this be the case? We expect to find a solution in which the photon is still approximately free, but 
has a mass gap. If this is the case, then entropically speaking the compact photon corresponds 
to a single free massive degree of freedom. But this possibility is contained in the expression for 
the entropy written in terms of the arbitrary kernels. Turning this argument upon its head, we 
can justify the conjecture a posteriori in the following fashion. Let us make the conjecture, carry 
out the minimisation and ask the question: Is the entropy thus obtained that of a single massive 
degree of freedom? As we shall see, it is. 

Putting everything together, the free energy in terms of G and h to order z is 


F = — 

4 J ( 2 t r ) 2 


^ r G -‘ ( i + 2A) +P 2 c(i-2Ar'-lAAA±^ 1! 

g z p z G z 



2 h 


'1- (1-4ft 2 ) 1 / 2 ' 

( i -( i -4 h 2 ) 1 ' 2 \y 

1 

(1 — 4 / t 2 ) 1 / 2 — (1 — 2h) 

log 

2h 

\ (1 — 4 / t 2 ) 1 / 2 — (1 — 2h) J J 


Minimising with respect to G, one obtains 

4n 2 7r 2 5z 


0 = -G _2 (l + 2 h) + p 2 ( 1 - 2/ij- 1 - 


r2 (1 + 2 h) 2 | 8n 2 TT 2 z(l + 2h) 2 


g 2 5G(p) 


d 2 k 


k 2 G 2 


g 2 p 2 G 3 


( 52 ) 


( 53 ) 


Minimising with respect to h, one obtains 


0 = G _1 +p 2 G(l — 2ft.) -2 


2 tt 2 7 t 2 Sz f 2 , (1 + 2 h) 2 8n 2 Tr 2 z(l + 2h) 

g 2 6h(p) J d k 2 G 2 + ^G 2 

2 T ri — (1 — 4 h 2 ) 1 / 2 ' 

+ (1 + 2h) 1 / 2 [\ — 2/i) 3 / 2 l0g [ 2 h _ 


From ( |33| ) one has the following relations 

Sz n 2 z 

Sh = ~2 g 2 p 2 G’ 
Sz n 2 z{ 1 + 2 h) 
SG = Ag 2 p 2 G 2 


( 54 ) 


( 55 ) 

( 56 ) 


We are interested in the weak coupling regime, where A /g 2 1, and the low temperature regime, 
where 0 < T < g 2 . If one evaluates the z terms in the minimisation equations (^3],^4|) in this regime, 
one finds that terms with h dominate over the others, which we discard. Thus, substituting 
for G from (^||) into ([hi]), one has 


1 - (1 - 4 / i 2 ) 1 / 2 
2 h 


(p 2 + m^/2 


= e 


( 57 ) 


which is the Boltzmann factor for a massive photon, where the mass is given by 

2 u 4 tt 2 z f (1 + 2 h ) 2 n 4 7r 3 A 2 2 

m '=— /"-w- = — 


(58) 


Here we we have used the approximations h ~ 0 and G _1 ~ k for large k. The form of the 
Boltzmann factor ensures that the entropy (which depends on h but not G) corresponds to the 
entropy of a single relativistic massive excitation of mass m. Thus our previous conjecture that 
compact 17(1) is (entropically speaking) a free theory is verified. 

®Because of the non-trivial gauge projection, the density matrix is no longer Gaussian. 
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It is now straightforward to solve for the kernels. One obtains 

2(p 2 +m 2 ) 1 / 2 

G - 1= ~ 


P 


1 + e 


(p 2 + m 2 ) 1 / 2 1 1 _ e _ a(p 3 +m=» ) i/» 
U, 2 -^ 2 ) 1 / 2 


if = 2 


(p 2 + m 2 ) 1 / 2 


2( p 2 +m 2)l/2 

1 - e t 


(59) 


In the limit as o 2 /A —> 0, m —> 0 and one obtains the expected kernels for the non-compact theory, 
cf. @ with w =p. 


5 Phase Structure of Compact U( 1) 

It is now straightforward to analyse the phase structure of the theory. From equations ( fbf ) and 
(p3|) the fugacity is given by 


A 2 n2lT 

* = A ' exp ~v. 

Now substitute for the kernels G and H from 

„ 2 _ r A 


pdp 


G~ L + H 

p2 


to obtain 


z — A 2 exp — - 


pdp 


1 


■ coth 


(p 2 + m 2 ) 1 / 2 


2 g 2 J Q ^ (p 2 + m 2 ) 1 /2 

The integral is straightforward. One obtains 

(p 2 + m 5 


2T 


(60) 


(61) 


A 2 '<■ " 

Z = A exp-— T 

2 r 


2 T log sinh ■ 


) 1/2 


A 


2 T 


- A 2 exp [A - 2T log m/T], 


= A 2 exp - 


25- 

n 2 nA 


m 

T 2 ) 


\ 2 s 2 


(62) 


Using the expression (| 
5 = z/ A 2 : 


with the two solutions 


for m yields a self-consistent equation for the dimensionless fugacity 


_nVA /n 4 7r 3 A 4 5\^ 
28 


5 = 0, 


(i / n 4 7r 3 A 4 \ 2 s 2 

S'- 2s 2 ’ = p 2 s 2 - 

1 p 4 t 2 ; 


(63) 


(64) 


The first solution is simple to comprehend. What about the second? The right hand side is, for 
0 < T < g 2 , a non-perturbatively small dimensionless quantity. At T = 0, one recovers the Kogan 
and Kovner result for the fugacity ||. As T increases from zero, the left hand side is z raised to 
a fractional power. This implies that the value of z which satisifes the equation is smaller and 
indeed goes continuously to zero as T goes to 2g 2 /n 2 n. More explicitly, writing t = tv"kT /2g 2 and 
taking the logarithm, one has that 


log 5 = — 


fit) 

1 -t' 


(65) 
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(66) 


where f(t) is a function of 0(A/g 2 ) which is both positive and smooth ( C °°) at t = 1. Thus, 

m 


z = exp — 


1 -t 


As t increases towards unity, 5 goes continuously to zero; above t = 1, z becomes large. 

What are we to make of all this? We have obtained two solutions which extremise the free 
energy. It is clear that the second solution minimises the free energy for temperatures below 
2g 2 /ro 2 7r (otherwise there would be no mass gap at zero temperature), but is unphysical (we only 
work to O(z) remember) for temperatures above 2g 2 /n 2 Tt and should be discarded. So at high 
temperatures we are left only with the solution z = 0. 

Hence there is a critical phase transition at 


T c = 2g 2 /n 2 n, 


( 67 ) 


from a phase with non-zero fugacity and mass gap (corresponding to an unbound plasma of 
monopoles and confinement of electric charges) to a phase with zero fugacity and mass gap (cor¬ 
responding to bound monopole-antimonopole pairs and deconfinement). 

Can we expand on the nature of this phase transition? Indeed we can. Consider again the low 
temperature solution (|66| ) for z. Not only does z go to zero as t tends to unity from below, but 
so also do all derivatives of z with respect to t, because / is smooth and the exponential falls off 
faster than any power of (1 — t). For example, 


.. dz rm 

inn — = hm e 1 ~ t 
t —i dt, t-> i 


' r 

i -t 


f ' 

(i-t) 2 _ 


= 0 , 


( 68 ) 


and so on. Thus we have that z and all its derivatives vanish on the low temperature side of the 
transition, as indeed they do on the high temperature side of transition (where z = 0). Thus z 
and all its derivatives are continuous at the transition. Moreover, this smoothness of z extends 
trivially to the free energy F, since any pathologies of F at t = 1 can only come from z. Thus, 
finally, we see that the free energy and all its derivatives are continuous at the transition: 


d m F 

dt m 


= 0. 


( 69 ) 


Thus, the transition is of infinite order according to the Ehrenfest classification. Since we do not 
know of another univerality class which has this property, we suggest that the phase transition 
does indeed belong to the BKT universality class, as expected. 


6 Correlation Functions and the Wilson Loop 


Having obtained the minimal variational kernels (|59|), it is now possible to calculate correlation 
functions in the compact 17(1) theory. For example, the magnetic field propagator is given by 
(following the arguments proceeding eq’n |5(]) 


( b(x)b(y )) = 


2trp' 


DADipDg 


A j ~ - g) 


tkidk 


A i ^ 7T d i(v-v) 

2 9 


-w 


( 70 ) 


where all quantities in the first set of parentheses take the argument x and so on. Hence, 


(b(x)b(y)) = idfcf (G- 1 - H)~\x - y), 

or in momentum space 


(b(p)b{-p)) 


( p 2 + TO 2 ) 1 / 2 

-coth 

2 


(p 2 + TO 2 ) 1 / 2 

2T 


( 71 ) 


( 72 ) 
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This is to be compared with the magnetic propagator for the non-compact theory at finite temper¬ 
ature calculated in the standard Lagrangian formalism. There (in Feynman gauge and neglecting 
decoupled ghosts) the Euclidean effective action is |lJ 



d 2 x -A fl ( d 2 + dl) Afj 


(73) 


obtained from the Minkowski action via the replacements t —> —it and Aq —> iA$. The magnetic 
propagator in terms of the Matsubara frequencies uj n = 2imT is then 


(■ B{x,T)B(y,r )) = E£L_. 


_fP_ P 2 i P ( X -y) 

(2 tt ) 2 P(p 2 +co 2 n ) e 


d 2 p 


p ,, p 
5 c «h-. 


(74) 


In the low temperature limit, the compact magnetic field propagator is simply (p 2 + 
and taking the Fourier transform (see e.g. 1141 ) one sees that the correlation functions 
e~m\x-y\ a |; l ar g e separations. In the high temperature limit, one has 


m 2 ) 1 / 2 /2 
decay as 


or in co-ordinate space 


(b(p)b(~p)) - T, 

(b(x)b(y)) —> TS(x — y). 


(75) 


(76) 


This too is to be expected: at high temperatures, fields undergo large thermal fluctuations but 
are uncorrelated at distinct space points. 

Strictly speaking, the correct propagator to consider is 


(b 2 (x)b 2 (y)) = 
1 


2trp' 


DADipDr) 


Cjj di 


A i - Yg d i^ ~ ^ 


i 2 r 


tkldk 


Al ~ 2 g dl( ' cp ~ ^ 


1 2 


W 


(77) 


since it is b 2 that is invariant under the full set of compact 17 ( 1 ) gauge transformations. There 
is a complication here in that there are terms containing 77 which appear to contribute to the 
propagator yet belong to an ill-defined path integral. The only 77 configurations which contribute 
to the path integral are those which are singular, since Cijdidjij = 0 when 77 is regular. But 77 only 
features in the path integral in the combination <p — 77 ; if we allow <p to be singular but insist that 
77 be regular, we generate all possible configurations p — rj (singular or otherwise) and at the same 
time ensure that the path integral is well defined. Via the Wick expansion, one obtains 


(b 2 (x)b 2 (y)} = 2 (b(x)b(y)} 2 , 


(78) 


which for small T falls off as e 2m ' [x v>i at large separations. 

Lastly, let us consider the spatial Wilson loop written in the form 

Was = ( e ilg fas A * dx *) = ( e il g!s BdS ) ; ( 79 ) 


where l takes integral values, ensuring that the Wilson loop is gauge invariant under compact 
transformations of the form B —> B + 2n/g. Then 


Wqs = j DAi exp —A(G 1 — H)A + ilg J tijdiAjcPx 


x / Dp exp 


1 ilf 

-—rd i p(G ~ 1 + H)dip+ — ® dipdxi . 

. 4 g 2 2 J dS J 


(80) 
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For a loop of large area, the first term is 


exp — J {B(x)B(y))d 2 xd 2 y 

= exp ^—S lim — H)~ x 

Z P —>o z 

l^Q f2 jyy 

= exp-Smooth—. (81) 

4 2T v ' 

The second term gives a vortex contribution of 0(z ) ~ 0(m 2 ) which is sub-dominant. So we have 
the apparently curious result that the area law behaviour persists for all T. Indeed, in the limit 
of small T, one has 


Was 


(82) 


which coincides with the result of Kogan and Kovner. In the limit of large T, one has 

J 2 2 m Q 

Was -> e-— (83) 

and so the string tension increases linearly with temperature. In fact, this too is to be expected. 
In general systems containing monopoles or vortices with some discrete magnetic symmetry group, 
the area law persists at high temperature, as discussed in jl5|] and confirmed by lattice simulations, 
e. g. @, 0 . Thus, the area law for the spatial Wilson loop is not a suitable confinement criterion 
at finite temperature . 

Remarkably, the same area law is found at finite temperature even for the non-compact theory. 
In this case the gauge-invariant Wilson loop can be written as 

W dS = (e i9 ^s A ^) = ( e i9l s BdSy (84) 

Proceeding as above and substituting for the non-compact kernels from ( |59| ) with m —> 0 yields 

Was-^e-^. (85) 

Thus, even for the Abelian U( 1) theory , large spatial Wilson loops show area law behaviour at all 
non-zero temperatures, with string tension proportional to the temperature. This result, which (to 
our knowledge) has not appeared in the literature thus far, is exact since the density matrix is 
exact in this case. It can be easily re-derived in the Lagrangian formalism using the propagator 

(ED- 


7 Discussion 

The variational analysis of the compact U(l) Hamiltonian at finite temperature has successfully 
reproduced all the results derived by dimensional reduction of the Lagrangian in |J. The expected 
critical transition has been identified (and at the expected temperature) and appears to belong 
to the expected universality class, viz. BKT. In addition, we have stressed the dependence of the 
critical temperature on the magnetic charge n of the fundamental monopoles in the system and 
shown that the area law is an inherent property of the Wilson loop at finite temperature. 

We identify two possible extensions. Firstly, it would be desirable to show unequivocally that 
the phase transition belongs to the BKT universality class (though we are not aware of any other 
phase transitions in two dimensional systems which are of infinite order). In order to do so, one 
needs to show that the couplings z and T obey renormalisation group equations which (after a 
suitable co-ordinate transformation in the parameter space of relevant couplings) have the BKT 
form. The equation (|b^) generates one RG equation, but it is not clear how to derive another 
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one in the Hamiltonian approach. Presumably, it would suffice to calculate some gauge-invariant 
Green’s function (and its scaling properties) in the theory using the minimal ansatz. 

Secondly, one could attempt a treatment of the full Georgi-Glasliow model using the variational 
approach. It has been shown jjj that in fact the deconfinement transition is driven by the unbinding 
of the Bosons at a temperature below that of the monopole binding. Could this be shown 
using a variational method? 

I would like to thank M. J. Bhaseen, I. I. Kogan and A. Kovner for their selfless assistance. 
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